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Abstract 

We investigate the confining properties of the QCD vacuum with Nf = 2 flavors 
of dynamical quarks, and compare the results with the properties of the quenched 
theory. We use non-perturbatively 0(a) improved Wilson fermions to keep cut-off 
effects small. We focus on color magnetic monopoles. Among the quantities we study 
are the monopole density and the monopole screening length, the static potential 
and the profile of the color electric flux tube. We furthermore derive the low-energy 
effective monopole action. Marked differences between the quenched and dynamical 
vacuum are found. 
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1 Introduction 

The dynamics of the QCD vacuum, and color confinement in particular, 
becomes more transparent and amenable to quantitative investigation in the 
maximally abelian gauge (MAG) [1,2]. In this gauge the relevant degrees of 
freedom are color electric charges, color magnetic monopoles, 'photons' and 
'gluons' [3]. The latter appear to become massive [4,5] due to a yet unresolved 
mechanism, resulting in an abelian effective theory at large distances. There 
is evidence that the monopoles condense in the low temperature phase of 
the theory [2,6], causing a dual Meissner effect which constricts the color 
electric fields into flux tubes, in accord with the dual superconductor picture 
of confinement. 

The dynamics of monopoles has been studied in detail in quenched lattice 
simulations. It turns out that in the MAG the string tension is accounted for 
almost entirely by the monopole part of the abelian projected gauge field [7,8], 
and that the low-energy effective monopole action is able to reproduce both the 
string tension and the low- lying glueball masses [9]. Furthermore, many of the 
non-perturbative features of the vacuum, such as the topological charge den- 
sity [10,11,12] and spontaneous chiral symmetry breaking [13], can be traced 
back to monopole excitations. 

Very little is known about the dynamics of monopoles in the full theory. So 
far the investigations have concentrated mainly on the static potential. While 
the effect of sea quarks is clearly visible at short distances, even for relatively 
heavy quark masses [14,15,16], no significant changes have been observed in 
the long-range behavior of the potential and the string tension. In contrast, 
the critical temperature of the chiral phase transition was found to depend 
noticeably on the mass of the dynamical quarks [17,18,19], which indicates 
that sea quarks have a visible effect on the non-perturbative properties of the 
vacuum as well. 

It will be interesting now to see how the microscopic properties of the 
vacuum react to the introduction of dynamical color electric charges. In this 
paper we shall study the effect of sea quarks on the dynamics of monopoles 
and the confining potential, and on the effective monopole action. The paper 
is organized as follows. In Section 2 we present the details of our simulations, 
as well as the gauge fixing procedure and abelian projection. In Section 3 
we discuss the gross properties of the vacuum, such as the monopole density 
and the magnetic screening length, and the static potential. Furthermore, the 
problem of Gribov copies is addressed. Section 4 is devoted to a detailed study 
of the static and dynamical properties of the color electric flux tube. In Section 
5 we derive the effective monopole action, employing an extended Swendsen 
method [20]. Finally, in Section 6 we conclude. Preliminary results of this work 
have been reported in Ref. [21]. 
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2 Simulation details 



Our studies are based on gauge field configurations with Nj = 2 flavors 
of dynamical quarks generated by the QCDSF-UKQCD collaboration, using 
the Wilson gauge field action and non-perturbatively 0(a) improved Wilson 
fermions [22]: 

S F = SP - )-ng c sw a 5 i ! {s)a- flu F IIU (s)ilj{s), (1) 

^ s 

where sP is the ordinary Wilson fermion action. Our data sample and run 
parameters are listed in Table 1. We will compare the results with the outcome 
of quenched simulations on lattices of similar size and lattice spacing. The 
parameters of our quenched runs are also given in Table 1. 

We fix the MAG [25] by maximizing the functional 

m = T^E(\Un(s,^\ 2 +\U 22 (s,^ + \U 33 (s,^) (2) 

S,fJ, 

with respect to local gauge transformations g of the lattice gauge field, 

U(s, /i) - U°(s, v) = g( S yU(s, fi)g(s + /}) . (3) 

To do so, we use a simulated annealing (SA) algorithm [26], in which the gauge 
transformed link variables U 9 are thermalized according to the probability 
distribution 

p(U°)=eMF[U 9 ]/T}, (4) 

where T is an auxiliary 'temperature' which is gradually decreased after every 
Monte Carlo sweep from T = 5 to T = 0.04. To do so, we use 7500 sweeps 
on the 16 3 32 lattice and 10000 sweeps on the 24 3 48 lattice. Every sweep 
consists of a heat bath update of each of the three SU(2) subgroups of the 
link matrices. After the final temperature has been reached, several local gauge 
transformations are applied until F[U] has attained its maximum value within 
machine precision. It is known [27,8] that the MAG is plagued by Gribov 
copies. This shows in the occurence of local maxima of F[U]. As a result, 
gauge non-invariant observables will, in general, depend on how close one gets 
to the global maximum. Although the SA algorithm performs much better 
than the iterative local maximization procedure used in SU (3) gauge theory 
so far, it is not always able to find the global maximum. We shall estimate the 
systematic error due to this problem in Section 3. 

The functional F[U] is invariant under local U(l) x U(l) gauge transfor- 
mations and global Weyl transformations. From the (gauge fixed) SU(3) link 
variables we derive abelian link variables [2,28] 

u(s,/jl) = diag(ui(s,f/),U2(s,Ai),u 3 (s,A0) » «»(«>A*) = exp(i^(s,//)) (5) 
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Table 1 

Parameter values of our dynamical (Nf = 2) [22] and quenched (Nf = 0) gauge 
field configurations. The improvement coefficient csw was computed in [23]. The 
quenched ro/a values have been taken from [24]. We have used ro = 0.5 fm to set 
the scale. 



with 



1 3 4 4 

6i(s, (i) = aig(Ua(s, n)) - - aig{Ujj(s, n)) ^ ^ G [--tt, -tt] . (6) 

j =1 



The abelian link variables u(s,fi) take values in U(l) x U(l). Under a general 
gauge transformation they transform as 

u(s, n) — > d(s)'u(s, /J,)d(s + jl) , 
d(s) = diag(exp(i«i(s)), exp(i oc 2 (s)), exp(-i(«i(s) + a 2 (s)))) • 

The monopole currents reside on links of the dual lattice and are denned 

by 

kiCs,ti) = ^- J2 axg( Wi (n))=0,±l,±2, (8) 

where Ui(n) is the product of abelian parallel transporters around the plaque- 
tte □, and f(s + jl, //) is the elementary cube perpendicular to the /i-direction 
with origin s + jl, with □ inheriting its orientation from df(s + //,//). (Note 
that we differ here from the original and correct normalization of the monopole 



4 



currents [3].) The phases are chosen such that 

5>rg(^(n)) = |arg(^(n)) -arg(^-(n))| < 2tt . (9) 

i 

Because of that, 

£>(* s ,a*) = o. (io) 

1=1 

3 Gross properties of the vacuum and static potential 

Let us first look at the global changes of the QCD vacuum upon introducing 
dynamical color charges. Our smallest quark masses are of the order of the 
strange quark mass. 

Monopole density 

The monopole currents, which are conserved, form clusters of closed loops 
on the dual lattice. In case of the pure SU (2) gauge theory it was observed 
that these clusters fall into two different classes [29,30]: 'small' (ultraviolet) 
clusters which are of limited extent in lattice units, and 'large' (infrared) 
clusters which percolate through the lattice and typically wrap around the 
boundaries. If the size of the lattice is large enough, a gap opens between 
the small and the large clusters, clearly separating the two. In general, each 
configuration accommodates at least one large cluster [30]. In Fig. 1 we show 
the histogram h(L) of monopole currents of length L on the 24 3 48 lattice. We 
observe two distinct clusters. For comparison we show the same quantity on 
the 16 3 32 lattice in Fig. 2. Both lattices have a similar lattice spacing and 
quark mass. On the smaller volume no gap is observed. For the long-distance 
properties of the vacuum all what matters is the existence of long, percolating 
monopole loops. Whether they combine to one cluster or not is of secondary 
importance. 

On smaller lattices we call a monopole cluster infrared if it forms the largest 
cluster or if the monopole loop wraps around the boundary. It has been shown 
in the pure SU(2) gauge theory [31] that the corresponding monopole density 
does not depend on the lattice volume and scales properly, while the total 
monopole density diverges in the continuum limit. 

We define the monopole density by 

^wEEN*^)!- (ii) 

In Fig. 3 we show the total monopole density as well as the density of monopoles 
belonging to the infrared cluster. In the former case the sum over s, fj, extends 
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Fig. 1. The histogram of closed monopole loops of length L in full QCD on the 
24 3 48 lattice at (3 = 5.29, k = 0.1355. The bin size is 200. 



2500 - 



25 - 



0.25 - 



0.025 




Fig. 2. The same as in Fig. 1, but on the 16 3 32 lattice at (3 = 5.29, k = 0.135. The 
bin size is 60. 



over the full lattice, while in the latter case it extends over the links of the 
infrared clusters only. The quenched result is entered at m^/rrip = 1. In the 
dynamical vacuum both densities are about a factor of two larger than in the 
quenched case. The total monopole density appears to increase with decreas- 
ing quark mass, while the density of infrared monopoles shows little variation, 
apart from the initial jump. 

How can one explain the increase of the monopole density in the dynam- 
ical vacuum? It has been known for some time that monopoles are induced 
by (anti-)instantons [10,11], at least partially. The fermion determinant in- 
troduces an attraction between instantons and anti-instantons, and the force 
increases with decreasing quark mass [32]. The effect is, very likely, that the 
vacuum becomes solidly packed with instantons and anti-instantons, while iso- 
lated instantons are suppressed. As a result the density of (anti-)instantons 
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Fig. 3. The monopole density in full and quenched QCD. The quenched results refer 
to P = 6.0 and are listed at m^jm, p = 1. 



increases, and consequently the density of monopoles. 

One should be aware that the total monopole density is not universal but 
depends, in general, on the action chosen, and that after integrating out the 
fermions the effective gluonic action may be more noisy than the Wilson gauge 
field action at the same lattice spacing. This would mainly affect the total 
monopole density, and to a lesser extent the infrared cluster. 



Static potential 

From the abelian projected link variables Ui(s,{i) we extract the abelian 
static quark-antiquark-potential. In order to improve the overlap with the 
ground state, we smear the space-like links according to 

Ui(s,j) -> aui(s,j) + ^Ui(s,k)ui(s + k,j)ui(s + ],k)^ . (12) 

We apply 30 smearing sweeps with a = 2. The abelian potential V a b(-R) is 
given by 

v (m v i ( <WAR,T)> \ .... 

where 

W ah (R,T) = l-ReTrWc, W c =\{u{s^), (14) 

and C is a (orientated) loop of spatial extent R and temporal extent T. 

The ratio on the r.h.s. of (13) reaches a plateau at T = 5, so that we will 
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Fig. 4. Comparison of the abelian and non-abelian potential in full QCD on the 
16 3 32 lattice at (3 = 5.29, k = 0.135. The solid lines are fits of the form (15). 
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Fig. 5. Same as in Fig. 5, but for the quenched theory at (5 = 6.0. 
take T = 5 throughout this section. We fit Kb(^) by the ansatz 



Kb(r) = K°b + ^ b r - 



«ab 



(15) 



The potential was calculated for on-axis and off-axis directions f = 1/ \[2 (1, 1, 0) 
and 1 / (1, 1, 1). At small r rotational symmetry is broken on the lattice, and 
we exclude the first four data points from our fits. The non-abelian static po- 
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Table 2 

The Coulomb term, the abelian and monopole part of the string tension, as well as 
the monopole screening length in full and quenched QCD. The quenched result is 
shown in the last row and refers to the 16 3 32 lattice at (3 = 6.0. 

tential was extracted from the non-abelian Wilson loop using a corresponding 
procedure. 

In Figs. 4 and 5 we show the abelian and non-abelian static potential for 
some data set. The self-energy contributions have been subtracted. The ratios 
of abelian and non-abelian string tensions are given in Table 2. As our deter- 
mination of the non-abelian string tension was not accurate enough, we took 
these numbers from the literature: \for§ = 1.142(5) in full QCD [33] (y/arQ 
depends only weakly on the dynamical quark mass) and y/aro = 1.16(1) in 
the quenched theory [34]. The abelian string tension turns out to be very close 
to the non-abelian one in full QCD, while in the quenched case it is noticeably 
smaller. But the ratio of <r a b to a may increase in the continuum limit [31]. 

The abelian link variables can be decomposed into a 'singular' monopole 
part and a photon part according to the definition [35,36]: 

6 l (s^)=6r n (s,ri+6l h (s^), (16) 
#rW) = 27Tj2D(s-s')V^m t (s',a,fi), (17) 

s' 

where D(s) = A _1 (s) is the lattice Coulomb propagator, is the lattice 
backward derivative, and rrii(s, /i, v) counts the number of Dirac strings pierc- 
ing the plaquette 

Ui(s, /I, v) = Ui(s, /i)ui(s + ft, v)u\{s + z>, fj,)u\(s, v) . (18) 

If one computes ki(*s,fi) from #™ on (s,/i) one recovers almost all monopole 
currents found by the definition (8), hence the notation monopole part [13]. 

Similarly, from Wilson loops composed of the monopole (photon) part of 
the link variables one can derive the monopole (photon) contribution to the 
static abelian potential. In Figs. 6 and 7 we show both contributions. The 
string tension of the monopole part a mon is given in Table 2, while the photon 
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Fig. 6. Decomposition of the abelian potential into monopole and photon parts on 
the 16 3 32 lattice at /3 = 5.29, k = 0.135 in full QCD. The solid lines are fits of the 
form (15). 




Fig. 7. The same as in Fig. 6, but for the quenched theory at (3 = 6.0. 



contribution to the potential is short-range. Within the uncertainties (of 0(a 2 ) 
corrections) all string tensions are very similar, in the dynamical as well as in 
the quenched theory. 

The next quantity we looked at is the magnetic screening length £. This is 
denned by the exponential decay of the magnetic flux <l>(r) through a sphere 
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Fig. 8. Data and fit of the magnetic flux on the 16 3 32 lattice in full and quenched 
QCD. 

of radius r around the monopole. On a periodic lattice this can be written 

$(r) = $ exp(- A) S mh(^— , (19) 

where L is the effective length of the box, which is taken to be a free parameter. 
Some numerical data are shown in Fig. 8, together with a fit of eq. (19). The 
fitted values of £ for all data sets are given in Table 2. The length L turns 
out to be slightly larger than the extent of the lattice, as expected. We notice 
that the screening length is about 30% lower in the dynamical vacuum as 
compared to the quenched case. This does not come unexpected. In a three- 
dimensional model of a monopole gas with screening [30] the abelian string 
tension turns out to be proportional to the product of monopole density and 
screening length, i.e. cr ab oc p£. Though this model is an oversimplification of 
the underlying dynamics, it is in qualitative agreement with our findings. 

Gribov copies 

We shall now try to quantify the error that is made by fixing to a local 
maximum of F[U] instead of the global one. We follow the procedure suggested 
in [8]. The test runs are done on the 16 3 32 lattice at (3 = 5.29, k = 0.135 using 
a total of O(50) configurations. We create 20 random gauge copies for each 
configuration, employing the SA algorithm. Then we randomly pick n gauge 
copies out of each set (of 20) and select the copy with the highest value of 
F[U] to compute our observable O. Obviously the result will depend on n, 
and the true result is obtained at n — > oo. The approach to n = oo may be 
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Fig. 9. The monopole density on the 16 3 32 lattice at (3 = 5.29, k = 0.135, and its 
dependence on the number of Gribov copies. 



fitted by [37] 

(0)(n) = <0>(°°) + — ■ (20) 

n 

In Fig. 9 we show the monopole density as a function of n. We see that p 
reaches a plateau at n ~ 10. By taking only one gauge copy into account one 
introduces a systematic error of the order of 3%. The effect of Gribov copies 
is slightly stronger in case of the abelian static potential, as can be inferred 
from Figs. 10 and 11. Here the systematic error is about 6%, while iterative 
gauge fixing might lead to a discrepancy of O(20%). 



4 Color electric flux tube 

Studies of the pure SU(2) gauge theory in the MAG [38,39,40] have shown 
that the expectation values of the static color electric field and the monopole 
currents satisfy, to a good accuracy, the classical equations of motion and 
dual Ampere's law, in agreement with the dual superconductor picture of 
confinement. 

In this section we present first results of the microscopic structure of the 
color electric flux tube in full QCD and in the pure 577(3) gauge theory. It 
is expected that long-range forces between quarks remain to exist in the dy- 
namical theory as well, because the color charge of quarks cannot be screened 
locally by Higgs scalars (made out of gluons). 
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Fig. 10. The abelian static potential on the 16 32 lattice at (3 = 5.29, re = 0.135, 
and its dependence on the number of Gribov copies. Also shown is the result of the 
iterative gauge fixing. 
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Fig. 11. The abelian string tension on the 16 3 32 lattice at (3 = 5.29, re = 0.135, 
and its dependence on the number of Gribov copies, together with the result of the 
iterative gauge fixing (o). 
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Observables 

We will primarily be concerned with local abelian operators 

O(s) = diag(CMs), 2 (s), O s (s)) E U(l) x 17(1) . (21) 

For C-parity even operators, such as the action density and the monopole 
density, the correlator of 0(s) with the abelian Wilson loop We is given by [41] 

For C-parity odd operators O, such as the color electric field and the monopole 
current, we have 

_ (Tr (O(s)Wc)) 

(G(s))w = — (^m) — ' (23) 

in analogy to the case of SU(2) [38,39,42,43]. 

The action density p™, the color electric field E™, the monopole current 
k w and the monopole density induced by the Wilson loop, are then given 
by 

Pa ( s ) = f S (diag(cos(#i(s, p, v)), cos(fl 2 (s, p, is)), cos(# 3 (s, p, v))))w , (24) 

where 9i(s,p, v) = arg(uj(s, p, v)) is the plaquette angle, 

Ef(s) = i (diag(^( S , 4, j), 9 2 (s, 4, j), 6 3 (s, 4, j))) w , (25) 



k w (*s, p) = 2vri (diag(A;i(*s, p), k 2 (*s, p), k 3 (*s, p))) w , (26) 

and 

P^(s)^lY,( d ^M*s,p)\,\U*s,p)\,\h(*s,p)\)) w , (27) 

respectively. Out of the three 'color' components of the observables only two 
are independent. In the following we shall take the average of the three com- 
ponents. 

As before, we take C to be a loop of spatial extent R and temporal extent 
T. The four corners of the loop are placed at (—i?/2,0,0,0), (-R/2, 0, 0, 0), 
(-R/2, 0, 0, T) and (R/2, 0, 0, T), and s 4 = 7/2 will be taken throughout this 
section. 

Abelian flux tube 

Let us first consider the profile of the abelian flux tube. We take R = 10, 
which on our lattices corresponds to a spatial separation of the static sources 
of ks 1 fm, and 7 = 6. We checked the 7 dependence of part of our results 
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Fig. 12. The action density p^(s)rQ of the Abelian flux tube as a function of 
x = si, y = S2 at z = S3 = on the 16 3 32 lattice in full (top) and quenched QCD 
(bottom) at (3 = 5.20, k = 0.1355 and (3 = 6.0, respectively. 

by comparing the numbers to T = 5 and found only insignificant changes, 
albeit for R = 6, which justifies our choice. The spatial links are smeared as 
described in (12). 

In the following we shall also use the notation x — s±, y — S2 and z = S3. In 
Figs. 12 and 13 we show the action density Pa(s) in full and quenched QCD, 
respectively. It appears that the action density in full QCD is higher than in 
the quenched case, while their shapes are quite similar. 

We estimate the width 5 of the abelian flux tube by fitting our data at 
x = to the function 

p A ( r± ) = const, exp ( - (r ± - ef/5 2 ) , (28) 

where r± is the distance of s from the line connecting the static sources at 
s 4 = T/2, and e is a displacement parameter of 0(a) accounting for a small 
shift of the true action density from its entry at s. We obtain 8 = 0.29(1) 
fm, both in full and quenched QCD. This is a surprisingly small number, 
much smaller than any hadron radius. It tells us, in particular, that already at 
interquark distances > 0.5 fm the string model of hadrons becomes effective. 
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In Figs. 14 and 15 we show the distribution of the color electric field in 
and around the flux tube. The electric field is purely longitudinal in a narrow 
region between the sources of diameter six lattice spacings and practically 
zero outside. We fit E™ at x, z = and for y/r > 0.5 to an exponential: 

= const. exp(-y/A) . (29) 

For the penetration length we find A = 0.15(1) fm in full QCD and A = 0.17(1) 
fm in the quenched case. Whether the flux tube indeed narrows down in full 
QCD needs to be confirmed with higher statistics. 

In Figs. 16 we show the monopole density Pm(s). We see again that out- 
side the flux tube the monopole density is more than two times larger in full 
QCD than in the quenched case. Inside the flux tube the monopole density 
is strongly suppressed. This indicates that the expectation value of the dual 
Higgs field vanishes inside the flux tube, in agreement with the dual super- 
conductor model of the vacuum. In this model we furthermore expect that 
the monopole currents forms a solenoidal (i.e. azimuthal) supercurrent which 
constricts the color electric field into flux tubes, thereby satisfying the dual 
Ampere law: 

k = V xE w . (30) 

In Fig. 17 we show the transverse components of the monopole current at 
x = in the (y,z) plane (i.e. perpendicular to the flux tube), and in Fig. 18 
we compare the l.h.s and r.h.s. of eq. (30). In the latter figure R = 6 was used in 
order to reduce the errors. We find that the dual Ampere law is approximately 
satisfied in both full and quenched QCD. So far this has only been verified in 
the pure SU(2) gauge theory [38,39]. 
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Fig. 14. Distribution of the color electric field in full (top) and quenched QCD 
(bottom) in the (x,y) plane for the same lattices as in Fig. 12. The magnitude of 
E w is indicated by the length of the arrows. 
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Fig. 16. The monopole density Pm(s)vq as a function of x, y at z = in full (top) 
and quenched QCD (bottom) for the same lattices as in Fig. 12. 




/ / . 

/ / / . 



Fig. 17. The solenoidal monopole current k w r^ in the (y,z) plane (i.e. perpendicular 
to the flux tube) at x = in full (left) and quenched QCD (right) for the same 
lattices as in Fig. 12. 
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Fig. 18. Test of dual Ampere's law in full (top) and quenched QCD (bottom) for 
the same lattices as in Fig. 12. 



Monopoles versus photons 

Further insight into the confinement mechanism can be obtained by probing 
the flux tube in terms of the monopole and photon part of the abelian gauge 
field separately. To do so, we simply have to replace the lattice abelian gauge 
field 9i(s, /i) in the various probes by 9™ on (s, jj) and Of h (s, //), respectively. We 
have done calculations in full QCD (here on the 24 3 48 lattice at (3 — 5.29, 
k = 0.1355) and in the quenched theory. To enhance the signal, the monopole 
and photon fields are smeared as before. Qualitatively, we find no difference 
between full QCD and the quenched theory. 

In Fig. 19 we show the action density of the monopole and photon part 
of the gauge field. We see that the action density originates almost entirely 



19 



Fig. 19. Monopole (top) and photon part (bottom) of the action density P^Tq on 
the (quenched) 16 3 32 lattice at (3 = 6.0. 




Fig. 20. Distribution of the monopole (top) and photon part of the color electric field 
E w (bottom) on the 16 3 32 lattice at (3 = 6.0. For better visibility the monopole 
part is enhanced by a factor of two relative to the photon part. 
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Fig. 21. Monopole and photon part of the color electric field in the (y,z) plane 
at x = (top) and parallel to the flux tube (bottom). 



from the monopole part, while the photon contributes a Coulomb field around 
the static charges only. The width of the flux tube is unchanged: 5 = 0.29(1) 
fm as before. In Figs. 20 and 21 we show the distribution of the color electric 
field. We see that the monopole part of the field has no sources. The sources 
show up in the photon part only. We furthermore see that outside the flux 
tube the monopole and photon parts of the electric field largely cancel, while 
they add in the interior of the tube. We have attempted to fit the photon field 
by a Coulomb ansatz. While the transverse component could be well fitted, 
we failed for the longitudinal component (i.e. parallel to the flux tube). 



String breaking 

In the presence of dynamical quarks we expect that the flux tube (string) 
breaks if the static charges are separated far enough. It has been estimated that 
this will happen at a distance of ~ 1.2 fm for our quark masses of m q 100 
Mev [14,44]. This does not mean that the string state vanishes from the spec- 
trum of the transfer matrix. It only ceases to be the state of lowest energy 
in the corresponding channel. In QCD string breaking has so far only been 
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Fig. 22. The monopole part of the color electric field E w on the 24 3 48 lattice. Only 
the region — 18 < x < 17, — 10 < y < 9 is shown. 

observed at finite temperature close to the deconfming phase transition [45], 
but never at zero temperature [14,15,16,46]. The recent finding [47] of string 
breaking from Wilson loops in the case of adjoint quarks in three-dimensional 
SU(2) gauge theory indicates that that this should be possible at zero tem- 
perature too. Though our prime motivation for the following investigation was 
to detect string breaking and shed some light on the dynamics that drives it, 
we like to stress that the unbroken string is of quite some interest as well from 
the point of view of the confinement problem. 

The calculations in this section are done on the 24 3 48 lattice at (3 = 5.29, 
k = 0.1355. On this lattice m n /m p « 0.7. The difference to our previous 
calculation is that we consider large Wilson loops. We take R = 18, which 
corresponds to a separation of the static charges of 1.6 fm, and T = 10. It 
is important to choose T large as well in order to increase the chance of string 
breaking [48]. In Fig. 22 we show the monopole part of the color electric field 
E w . The restriction to the monopole part allowed us to obtain a clean signal 
even at x — 0. We do not observe any sign of string breaking. Furthermore, 
the flux tube does not show any broadening effect. A fit of the form (28) gives 
5 = 0.30(3) fm at x = 0. The same value was found for R = 10. If the long- 
range properties of QCD were described by the Nambu-Goto effective string 
action, we would have expected the transverse extension of the flux tube to 
increase logarithmically with R [49]. A similar observation has been made in 
the pure SU(2) gauge theory [50]. Perhaps the abelian flux tube is described 
by the Ramond effective string action instead, which does not give rise to any 
broadening effect [51]. 

It is perhaps not surprising that we do not observe any string break- 
ing (yet) [52]. String breaking is expected to occur if exp(— 2E s i(R + T)) > 
exp(— crRT), where E s i is the binding energy of the static- light meson. This is 
only the case if T > 3 fm [16]. 
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5 Effective monopole action 



We have seen that the vacuum undergoes several changes if dynamical color 
electric charges are introduced. We shall study now how this will affect the 
effective monopole action. 

There are three types of monopole currents (8), of which two are indepen- 
dent. For simplicity we take into account only one of them, thus integrating 
out the other two [53] . For the time being, we assume the form of the effective 
monopole action in full QCD to be the same as in the quenched theory [54]. It 
is composed of 27 types of two-point interactions, one four-point interaction 
and one six-point interaction: 

29 

SW^GiSiik), (31) 

i=l 

where are the coupling constants which need to be determined. This we 
will do by employing an extended Swendsen method [20]. In this section we 
shall write k(s,/i) instead of k(*s,/i) for the sake of simplicity. Explicitly, we 
then have: 

Two-point interaction for parallel currents 

S l (k) = ^2T, k (s,fi)k l (s,fi), * = 1, •••,25, (32) 

S [M=l 

where the ki(s,/i) are given in Table 3. 
Two-point interaction for orthogonal currents 

S 2 %{k) = ^^k(s^)k(s - [i-2v,v) , (33) 
S2i{k) = Y. E k{s,n)k{s-jl-2v-2a,a). (34) 

s fi^u^cr 

Four-point interaction 

S28(k) = Y.(i: k(s,n) 2 ) • (35) 

s \ M =-4 / 

Six-point interaction 

S 2 9{k) = E ( E Hs,u) 2 ) • (36) 

s j 

The calculations are done on the 24 3 48 lattices listed in Table 1, both in 
full and quenched QCD. After fixing the gauge fields to the MAG, we employ 
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Table 3 

The monopole currents ki(s,p) entering eq. (32). 

a type-II block spin transformation [55] with up to n = 4 blocking steps. The 
final outcome is an action at the physical length scale b = na. We believe that 
the monopole action is effective at scales 0.4 < b < 0.8 fm. 

In Fig. 23 we show the self-coupling G\ as a function of b. We see that 
in full QCD Gi is systematically smaller than in the quenched theory for all 
values of b. In Fig. 24 we plot the self-coupling G\ as a function of m^/rrip for 
our smallest b value. We find that G\ decreases with decreasing quark mass. 

A necessary condition for monopole condensation is G\ < In 7. This is 
achieved for b > 0.27 fm in full QCD and for b > 0.35 fm in the quenched 
theory. In Fig. 25 we show the coupling constants G 2 and G 3 . We see that 
in full QCD G 2 is systematically smaller than in the quenched theory for 
all values of b, as in the case of G±, while G3 shows the opposite behavior. 
The other coupling constants show little difference between full and quenched 
QCD. In Fig. 26 we show, as an example, G28 and G29. Because the magnetic 
charge of the monopoles is in almost all cases ±1, and Gi is the dominant 
coupling constant, the monopole action can be approximated by Gi L, where 
L is the length of the monopole loop. As a result, a smaller self-coupling 
G\ will give rise to a larger value of L and a larger monopole density, and 
vice versa. This is to say that both observations, namely that the monopole 
density increases and the self-coupling decreases with decreasing quark mass, 
are consistent with each other. 
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Fig. 23. The monopole self-coupling G±, as a function of the physical length scale 
b. The symbols are: n = 1 (•), n = 2 (■), n = 3 (♦), n = 4 (a) in full QCD, and 
n = 1 (O), n = 2 (□), n = 3 (0), n = 4 (A) in quenched QCD. 
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Fig. 24. The dependence of G\ on the ratio m n /m p for b = 0.09 fm. The entry at 
m n /m p = 1 corresponds to the quenched theory. 
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Fig. 25. The coupling constants G2 and G3 as a function of b. The symbols are as 
in Fig. 23. 



26 



00 






0.001 



0.0005 - 



ON 





-0.0005 




Fig. 26. The coupling constants G28 and G29 of the four-point and six-point inter- 
actions, respectively, as a function of b. The symbols are as in Fig. 23. 
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Fig. 27. The two-point coupling G of the monopole currents k(s,fi) and k(s',/j,) 
as a function of distance d = \ /Yla( s ^ ~ s p 2 m ^ un (*) an d quenched QCD (■), 



compared to the Coulomb propagator (<>)• 
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Fig. 28. The self-coupling G se if and the Coulomb coupling g m as a function of b. 
The symbols are as in Fig. 23. 
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To shed further light on the dynamics of the monopoles we have looked at 
the coupling G of two (n = 1) parallel monopole currents, k(s, //) and k(s', //), 
as a function of the distance d = ^Jj2n( s ^ ~ s '^) 2 between them. In Fig. 27 
we show G together with the lattice Coulomb propagator. We see that at 
distances d > 2 the interaction becomes weaker than Coulomb in both full 
and quenched QCD. This is consistent with the screening effect discussed in 
Section 3. We do not see any difference between full QCD and quenched QCD 
though. 

We may parameterize the effective monopole action by 

= ^Coulomb + ^self + S4_ po int + S^-point + S'add, (37) 

where 

2 4 

Scoulomb = yEE k ( S , ») A" 1 (S ~ S')k(s', /i) , 
Z s,s> n=l 



S'self — G se \f Si(k) , 
S4-point — C4_ po i nt S28(k) , (38) 
S 6 -point — G^-point S^fc) , 

and S a dd includes 12 additional two-point interaction terms. In Fig. 28 we show 
the self-coupling G S eif and the Coulomb coupling g m . It is interesting to see 
that in full QCD G sc \{ is smaller than in the quenched case, while the Coulomb 
coupling is almost unchanged. Corrections to the Coulomb interaction are 
found to be very small in the infrared region. 



6 Conclusions 

We have performed a detailed study of the dynamics of the QCD vacuum, 
thereby focussing on the abelian degrees of freedom in the MAG. Our main 
objective was to find out how the vacuum reacts to the introduction of dy- 
namical color electric charges (quarks). The monopole density was found to 
increase by more than a factor of two if we decrease the quark mass from 
m n /m p = 1 (the quenched limit) to m n /m p ps 0.6, both for the total number 
of monopoles and for the monopoles in the infrared clusters. Related to that, 
we found that the magnetic screening length decreased by 30% over that range. 
The string tension, the static potential and the structure of the flux tube, on 
the other hand, remained almost the same. We verified the dual Ampere law 
in full QCD and in the pure SU(3) gauge theory. This result lends further 
support to the dual superconductor model of the vacuum in full and quenched 
QCD. The width of the abelian flux tube was found to be 5 = 0.29(1) fm in 
both cases. Another characteristic feature of the flux tube is the penetration 
length. We obtained A = 0.15(1) fm in full QCD and A = 0.17(1) fm in the 
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quenched case. This results in a dual photon mass of 1.3(1) GeV and 1.2(1) 
GeV, respectively. Decomposing the abelian gauge field into monopole and 
photon parts allowed us to study flux tubes up to a length of ~ 1.6 fm in full 
QCD. No signal of string breaking was found. Comparing flux tubes of various 
lengths R, it turned out that the width of the flux tube does not depend on 
R, contrary to the prediction of the Nambu-Goto effective string theory. The 
effective monopole action was determined. In full QCD the monopole self- 
coupling was found to be systematically smaller than in the quenched theory. 
The main contributions to the effective monopole action are found to be the 
self-interaction and the Coulomb interaction. 
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